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' Abstract 

CN ■ 

The Eyring-Kramers law describes the mean transition time of an overdamped Brown- 
^ . ian particle between local minima in a potential landscape. In the weak-noise limit, the 

' transition time is to leading order exponential in the potential difference to overcome. 

This exponential is corrected by a prefactor which depends on the principal curvatures 
of the potential at the starting minimum and at the highest saddle crossed by an opti- 
mal transition path. The Eyring-Kramers law, however, does not hold whenever one 
or more of these principal curvatures vanishes, since it would predict a vanishing or 
infinite transition time. We derive the correct prefactor up to multiplicative errors 

that tend to one in the zero-noise limit. As an illustration, we discuss the case of a 
• I ' 

. symmetric pitchfork bifurcation, in which the prefactor can be expressed in terms of 

' modified Bessel functions, as well as bifurcations with two vanishing eigenvalues. The 

corresponding transition times are studied in a full neighbourhood of the bifurcation 
point. These results extend work by Bovier, Eckhoff, Gayrard and Klein |BEGK 04J. 
who rigorously analysed the case of quadratic saddles, using methods from potential 
' theory. 
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1 Introduction 



Consider the stochastic differential equation 

dxt = -VV{xt) dt + V2e dWt , (1.1) 



where y : M"' — > R is a confining potential. The Eyring-Kramers law ( |Eyr35 IKra40] ) 



describes the expected transition time r between potential minima in the small-noise limit 
e — > 0. In the one-dimensional case (d = 1), it has the following form. Assume x and y 
are quadratic local minima of V, separated by a unique quadratic local maximum z. Then 
the expected transition time from x to y satisfies 

E-|r| ~ ^[Vizyvix)]/e ^ (1^2) 

^ ^ VF"(x)|y"(z)| 

In the multidimensional case (d ^ 2), assume the local minima are separated by a unique 
saddle z, which is such that the Hessian \/'^V{z) admits a single negative eigenvalue Ai(2;), 
while all other eigenvalues are strictly positive. Then the analogue of (|1.2p reads 



E-/rl ~ 27r / det(W(z)) [v(z)-V(x)]/e ^ 3^ 

^ |Ai(z)|\/det(W(x))'' ■ ^^-"^^ 

This expression has been generalised to situations where there are several alternative 
saddles allowing to go from x to y, and to potentials with more than two minima. 

A long time has elapsed between the first presentation of the formula (jl.Sp by Eyring 
|Eyr35| and Kramers |Kra40j and its rigorous mathematical proof (including a precise 
definition of what the symbol "~" in (jl.3p actually means). The exponential asymp- 
totics were proved to be correct by Wentzell and Freidlin in the early Seventies, using 
the theory of large deviations |VF69l IVFTOi IFW98] . While being very flexible, and allow- 
ing to study more general than gradient systems like (jl.ip . large deviations do not allow 
to obtain the prefactor of the transition time. An alternative approach is based on the 
fact that mean transition times obey certain elliptic partial differential equations, whose 
solutions can be approximated by WKB-theory (for a recent survey of these methods, 
see [KolOOj ). This approach provides formal asymptotic series expansions in e, whose 
justification is, however, a difficult problem of analysis. A framework for such a rigor- 
ous justification is provided by microlocal analysis, which was primarily developed by 
Helffer and Sjostrand to solve quantum mechanical tunnelling problems in the semiclas- 
sical limit |HS84t IHS85bl IHSSSaj IHS85cj . Unfortunately, it turns out that when trans- 
lated into terms of semiclassical analysis, the problem of proving the Eyring-Kramers 
formula becomes a particularly intricate one, known as "tunnelling through non-resonant 
wells" . The first mathematically rigorous proof of (11. 3p in arbitrary dimension (and its 
generalisations to more than two wells) was obtained by Bovier, Eckhoff, Gayrard and 
Klein [BEGK04] . using a different approach based on potential theory and a variational 
principle. In [BEGK04] . the Eyring-Kramers law is shown to hold with a ~ 6 meaning 
a = 6(1 + ©(e^/^llogel)). Finally, a full asymptotic expansion of the prefactor in powers 
of e was proved to hold in |HKN04| IHNOSj . using again analytical methods. 

In this work, we are concerned with the case where the determinant of one of the 
Hessian matrices vanishes. In such a case, the expression (II. 3p either diverges or goes 
to zero, which is obviously absurd. It seems reasonable (as has been pointed out, e.g., 
in |Ste05j ) that one has to take into account higher-order terms of the Taylor expansion 
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of the potential at the stationary points when estimating the transition time. Of course, 
cases with degenerate Hessian are in a sense not generic, so why should we care about 
this situation at all? The answer is that as soon as the potential depends on a parameter, 
degenerate stationary points are bound to occur, most noteably at bifurcation points, i.e., 
where the number of saddles varies as the parameter changes. See, for instance, [BFG07al 
IBFGOTb] for an analysis of a naturally arising parameter-dependent system displaying a 
series of symmetry-breaking bifurcations. For this particular system, an analysis of the 
sub exponential asymptotics of metastable transition times in the synchronisation regime 
has been given recently in [BBMOO j. with a careful control of the dimension-dependence 
of the error terms. 

In order to study sharp asymptotics of expected transition times, we rely on the 
potential-theoretic approach developed in |BEGK04l IBGK05] . In particular, the expected 
transition time can be expressed in terms of so-called Newtonian capacities between sets, 
which can in turn be estimated by a variational principle involving Dirichlet forms. The 
main new aspect of the present work is that we estimate capacities in cases involving 
nonquadratic saddles. 

In the non-degenerate case, saddles are easy to define: they are stationary points at 
which the Hessian has exactly one strictly negative eigenvalue, all other eigenvalues being 
strictly positive. When the determinant of the Hessian vanishes, the situation is not so 
simple, since the nature of the stationary point depends on higher-order terms in the Taylor 
expansion. We thus start, in Section [2l by defining and classifying saddles in degenerate 
cases. In Section [3l we estimate capacities for the most generic cases, which then allows 
us to derive expressions for the expected transition times. In Section HI we extend these 
results to a number of bifurcation scenarios arising in typical applications, that is, we 
consider parameter-dependent potentials for parameter values in a full neighbourhood of 
a critical parameter value yielding non-quadratic saddles. Section [5] contains the proofs of 
the main results. 

Acknowledgements: We would like to thank Bastien Fernandez for helpful discus- 
sions and Anton Bovier for providing a preliminary version of |BBM09) . BG thanks 
the MAPMO, Orleans, and NB the CRC 701 Spectral Structures and Topological Methods 
in Mathematics at the University of Bielefeld, for kind hospitality. Financial support by 
the French Ministry of Research, by way of the Action Concertee Incitative (ACI) Je- 
unes Chercheurs, Modelisation stochastique de systemes hors equilibre, and the German 
Research Council (DFG), by way of the CRC 701 Spectral Structures and Topological 
Methods in Mathematics, is gratefully acknowledged. 

2 Classification of nonquadratic saddles 

We consider a continuous, confining potential F : M'^ ^ M, bounded below by some 
ao € K and having exponentially tight level sets, that is, 

f e-^(^)/^ dx ^ C7(a) e"'*/" Va > oq , (2.1) 

with C{a) bounded above and uniform in e ^ 1. We start by giving a topological definition 
of saddles, before classifying saddles for sufficiently differentiable potentials V. 
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Figure 1. Examples of potentials and gates, (a) G{A, B) ~ {z}. (b) G{A, B) — {zi, Z2}. 
(c) G{A, B) = {zi} or {2:2}- Here curves show level lines, shaded areas indicate potential 
wells and the star marks a potential maximum. 



2.1 Topological definition of saddles 

We start by introducing the notion of a gate between two sets A and B. Roughly speaking, 
a gate is a set that cannot be avoided by those paths going from A to B which stay as low 
as possible in the potential landscape. Saddles will then be defined as particular points 
in gates. 

It is useful to introduce some terminology and notations: 

• For x,y G M"^, we denote by 7 : x ^ y a path from x to y, that is, a continuous 
function 7: [0, 1] ^ M*^ such that 7(0) = x and 7(1) = y. 

• The communication height between x and y is the highest potential value no path 
leading from x to y can avoid reaching, even when staying as low as possible, i.e., 

V{x,y)= inf sup V{j{t)) . (2.2) 

Note that V{x,y) ^ V{x) V V{y), with equality holding, for instance, in cases where 
X and y are "on the same side of a mountain slope" . 

• The communication height between two sets A, S C M is given by 

V{A,B)= inf V{x,y). (2.3) 

We denote by g{A,B) = {z G M'^: V{z) = V{A,B)} the level set oiV{A,B). 

• The set of minimal paths from A to S is 

ViA,B) = \^:x^y x e A,y e B, sup V{-fit)) = V{A, B)} . (2.4) 

te[o,i] J 

The following definition is taken from [BEGKOl] . 

Definition 2.1. A gate G{A,B) is a minimal subset of g{A, B) such that all minimal 
paths 7 G V{A,B) must intersect G{A,B). 

Let us consider some examples in dimension d = 2 (Figure [T|) : 

• In uninteresting cases, e.g. for A and B on the same side of a slope, the gate G{A, B) 
is a subset of ^ U i?. We will not be concerned with such cases. 

• If on the way from ^ to -B, one has to cross one "mountain pass" z which is higher 
than all other passes, then G{A,B) = {z} (Figure [TK)- 
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Figure 2. Examples of potentials and valleys. In case (b), a; is a saddle. 



• If there are several passes at communication height V{A, B) between A and i?, between 
which one can choose, then the gate G{A, B) is the union of these passes (Figure [Da) . 

• If when going from A to B, one has to cross several passes in a row, all at communica- 
tion height V{A, B), then the gate G{A, B) is not uniquely defined: any of the passes 
will form a gate (Figure [it) . 

• If A and B are separated by a ridge of constant altitude V{A, B), then the whole ridge 
is the gate G{A,B). 

• If the potential contains a flat part separating A from B, at height V{A, B), then any 
curve in this part separating the two sets is a gate. 

We now proceed to defining saddles as particular cases of isolated points in gates. 
However, the definition should be independent of the choice of sets A and B. In order to 
do this, we start by introducing notions of valleys (cf. Figure [2]): 

• The closed valley of a point x € M is the set 

CV{x) = {y eR'^:V{y,x) = V{x)} . (2.5) 

It is straightforward to check that CV{x) is closed and path-connected. 

• The open valley of a point x G M is the set 

OV{x) = {ye CV{x) : V{y) < V{x)} . (2.6) 

It is again easy to check that OV{x) is open. Note however that if the potential 
contains horizontal parts, then CV{x) need not be the closure of OV{x) (Figure [It). 
Also note that OV{x) need not be path-connected (Figure [2)3) . We will use this fact 
to define a saddle. 

Let Bs{x) = {y G M*^: \\y — x\\2 < e} denote the open ball of radius e, centred in x. 

Definition 2.2. A saddle is a point z G M"^ such that there exists e > for which 

1. OV{z) n B£{z) is non-empty and not path-connected. 

2. {pV{z)\J {z\)r\Bs{z) is path- connected. 

The link between saddles and gates is made clear by the following two results. 

Proposition 2.3. Let z he a saddle. Assume OV{z) is not path- connected]^ and let A and 
B belong to different path- connected components ofOV{z) nBe{z). Then z £ G{A,B). 

^We need to make this assumption globally, in order to rule out situations where z is not the lowest 
saddle between two domains. 
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Proof: Choose points x ^ A and y Si B and a path 7: A ^ B. Since A and B belong to 
different path-connected components of OV{z), the path 7 must leave OV(z), which implies 
supjgjo,!] ^(7(*)) ^ ^i^)- Since {OV{z)U{z})r]B^{z) is path-connected, we can find a path 
'y: A^ B staying all the time in this set, and thus for this path, sup^gjo,!] ^(7(0) — ^(-^)- 
As a consequence, the communication height V{x,y) equals V{z), i.e., 7 belongs to the 
set V{A,B) of minimal paths. Since OV{z) is not path-connected, we have found a path 
7 G V{A, B) which must contain z, and thus z G G{A, B). □ 

Proposition 2.4. Let A and B he two disjoint sets, and let z G G{A,B). Assume that 
z is isolated in the sense that there exists e > such that B*{z) :=Be{z) \ {z} is disjoint 
from the union of all gates G{A,B) between A and B. Then z is a saddle. 

Proof: Consider the set 

^= U U ^(*) (2-7) 
jeV{A,B) te[o,i] 

of all points contained in minimal paths from A to B. We claim that D = CV{z). 

On one hand, if a; G -D then there exists a minimal path from AtoB containing x. We 
follow this path backwards from x to A. Then there is a (possibly different) minimal path 
leading from the first path's endpoint in A through z to B. By gluing together these paths, 
wc obtain a minimal path connecting x and z. This path never exceeds the potential value 
V(z), which proves V{x, z) = V{z). Thus x G CV{z), and D C CV{z) follows. 

On the other hand, pick y G CV{z). Then there is a path ji: y —>■ z along which 
the potential does not exceed V{z). Inserting this path (twice, going back and forth) in 
a minimal path 7 G V{A,B) containing z, we get another minimal path from A to B, 
containing y. This proves y E D, and thus the inverse inclusion CV{z) C D. 

Now pick X £ A and y £ B. There must exist a minimal path 7: x ^ y, containing z, 
with the property that V is strictly smaller than V{z) on 7([0, l])nB*{z), since otherwise 
we would contradict the assumption that z be isolated. We can thus pick x' on 7 between 
X and z and y' on 7 between z and y such that V{x') < V{z) and V{y') < V{z). Hence 
we have x',y' G OV{z) and any minimal path from x' to y' staying in Be{z) has to cross 
z ^ OV{z). This shows that OV{z)nBe{z) is non-empty and not path-connected. Finally, 
take any x,y E OV{z)nB^{z) C D. Then we can connect them by a path j 3 z, and making 
£ small enough we may assume that V is strictly smaller than V{z) on 7([0, 1]) H B*{z), 
i.e., 7([0, 1]) \ {z} C OV{z). This proves that (C'V(z) U {z})nBe{z) is path-connected. □ 

2.2 Classification of saddles for differentiable potentials 

Let us show that for sufficiently smooth potentials, our definition of saddles is consis- 
tent with the usual definition of nondegenerate saddles. Then we will start classifying 
degenerate saddles. 

Proposition 2.5. Let V be of class , and let z he a saddle. Then z is a stationary 
point ofV, i.e., W{z) = 0. 

Proof: Suppose, to the contrary, that 'SJV{z) 7^ 0. We may assume z = and V{z) = 0. 
Choose local coordinates in which Vy(0) = (a, 0, . . . , 0) with a > 0. By the implicit- 
function theorem, there exists a differentiable function h : R'^"^ M*^ such that all 
solutions of the equation V{x) = in a small ball ^^(O) are of the form xi = h{x2, ■ ■ ■ , x^). 
Furthermore, V{e, 0, . . . , 0) = ae -|- 0(e) is positive for e > and negative for e < 0. By 
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continuity, V{x) is positive for xi > h{x2, ■ ■ ■ , x^) and negative for xi < h{x2, ■ ■ ■ jXd), 
showing that OV{z) D -^^(O) is path-connected. Hence z is not a saddle. □ 

Proposition 2.6. Assume V is of class C^, and let z be a saddle. Then 

1. the Hessian V'^V{z) has at least one eigenvalue smaller or equal than 0; 

2. the Hessian V'^V{z) has at most one eigenvalue strictly smaller than 0. 

Proof: Denote the eigenvalues of V^V{z) hy Xi ^ X2 ^ ■ ■ ■ ^ Xd- We may again assume 
that z = and V{0) = 0, and choose a basis in which the Hessian is diagonal. Then 

1 

V{x) = -^X,xj + oi\\x\\l). (2.8) 

i=l 

1. Assume, to the contrary that Ai > 0. Then V > near z = 0, so that OV{z) = 0, and 
2; = is not a saddle. 

2. Suppose, to the contrary, that Ai ^ A2 < 0, and fix a small S > 0. Since 

V{x) = ~\X,\xl - l|A2|xi + ^^^i + <^(ll^lli) ' (2-9) 

1=3 

we can find an e = e(5) G (0, S) such that for any fixed (xs, . . . , x^) of length less than 
e, the set {{xi, X2) : + X2 < (5^, V{x) < 0} is path-connected (topologically, it is an 
annulus). This implies that {(xi, . . . , Xd) : xf -|- x^ < (5^ , V{x) < 0, ||(x3, . . . , Xd)\\2 < 
is also path-connected. Hence OV{z) fl Bs{z) is path-connected. □ 

Proposition 2.7. Assume V is of class C^, and let z be a nondegenerate stationary point, 
i.e. such that det(V^y(2;)) 7^ 0. Then z is a saddle if and only ifV'^V{z) has exactly one 
strictly negative eigenvalue. 

Proof: Denote the eigenvalues of V^F(z) hy Xi ^ X2 ^ ■ ■ ■ ^ Xd. If z is a saddle, then 
the previous result implies that Ai < < A2. Conversely, if Ai < < A2, we have 

Vix) = -^|Ai|x? + i ^ A,xf + o(||x||^) . (2.10) 

i=2 

Thus for fixed small (x2, . . . ,Xd), the set {xi : |xi| < e, V{x) < 0} is not path-connected, 
as it does not contain 0. Thus OV{z) Ci I3£{z) is not path-connected (it is topologically 
the interior of a double cone). However, for xi = 0, adding the origin makes the set 
path-connected, so that {OV{z) fl Be{z)) U {z} is path-connected. □ 

We can now classify all candidates for saddles in the following way. Let Ai ^ A2 ^ 
■ ■ ■ ^ Xd he the eigenvalues of the Hessian V'^V{z) of a stationary point z, arranged in 
increasing order. Then the following cases may occur: 

1. Ai < 0: 

(a) A2 > 0: z is a nondegenerate saddle. 

(b) A2 = 0: 

i. A3 > 0: z is a singularity of codimension 1. 

ii. A3 = 0: 
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A. A4 > 0: 2; is a singularity of co dimension 2. 

B. A4 = 0: z is a singularity of codimension larger than 2. 

2. Ai = 0: 

(a) A2 > 0: z is a singularity of codimension 1. 

(b) A2 = 0: 

i. A3 > 0: 2; is a singularity of codimension 2. 

ii. A3 = 0: z is a singularity of codimension larger than 2. 

One can of course push further the classification, including all singularities up to 
codimension d. 

2.3 Singularities of codimension 1 

We assume in this subsection that the potential V is of class and that z is a stationary 
point of V with the Hessian \7'^V{z) having one vanishing eigenvalue. We may assume z = 
and V{z) = 0. According to Proposition 12.61 there are two cases to be considered: 

1. Ai < 0, A2 = and < A3 ^ • • • ^ A^. 

2. Ai = and < A2 ^ • • • ^ Arf. 

It will be convenient for the purpose of this subsection to relabel the first two eigenvalues 
in such a way that Ai = 0, while A2 7^ can be positive or negative and to choose a 
basis in which V'^V{z) = diag(0, A2, • • • , A^^). For potentials V of class and ii, . . . , v G 
{1, . . . , d}, we introduce the notation 

Vt,..,^ = — . Viz), nj = #{k:ik=j} , (2.11) 

Til. ■ ■ ■ IT'd' ClXi^ OXi^ 

with the convention that the z^'s are always in increasing order. 

Under the above assumptions, the potential admits a Taylor expansion of the form 

1 

= i^'^^ixi + ^ VijkXiXjXk+ ^ VijkiXiXjXkXi + o{\\x\\2) , {2.12) 

1=2 l^i^j^k^d l^i^j^k^l^d 

and the theory of normal forms allows to simplify this expression. Indeed, there exists a 
change of variables x = y + g{y), with g a polynomial function, such that the potential 
expressed in the new variables has as few as possible terms of low order in its Taylor 
expansion. In general, only a few so-called resonant terms cannot be eliminated, and are 
thus essential for describing the local dynamics. 

Proposition 2.8. There exists a polynomial change of variables x = y + g{y), where g is 
a polynomial with terms of degree 2 and 3, such that 

1 

V{y + g{y)) = - ^ A.yf + C^yl + C^y\ + o(||y||^) , (2.13) 

i=2 

where 

C3 = Fill , C4 = ^1111 - E 4^ • (2-14) 

2, =2 

The proof uses standard normal form theory and will be given in Appendix [Aj Note 
that if V is of class C^, then (|2.13p holds with odlylH) replaced by 0(112/112)- Let us now 
apply the result to derive an easy to verify necessary condition for a point z; to be a saddle. 
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Corollary 2.9. 

1. Assume A2 < 0. Then the point z is 

• a saddle if C3 = and > 0; 

• not a saddle if C3 ^ or C4 < 0. 

2. Assume A2 > 0. Then the point z is 

• a saddle i/ C3 = and < 0; 

• not a saddle if C3 ^ or C4 > 0. 

Proof: Consider first the case C3 7^ 0. For simplicity, let us restrict to d = 2. In a 
neighbourhood of z = 0, any solution to the equation V{y) = must satisfy 

yl = -^-^y!-^-^yf + oi\\y\\i). (2.15) 

Thus, solutions exist for yi with yiC^/\2 < 0. Plugging the ansatz 



y2 = ±^-^yl[l + r2{yi)] (2.16) 

into the relation V{y) = 0, dividing by yf and applying the implicit-function theorem to 
the pair (r2 ,yi) in the resulting equation shows that there is a unique curve through the 
origin on which the potential vanishes. Since for y2 = 0, the potential has the same sign as 
yi, we conclude that is not a saddle. Now just note that the proof is similar in dimension 
d>2. 

Consider next the case C3 = and A2C4 > 0. If A2 > 0, one sees that the origin 
is a local minimum. If A2 < 0, then for fixed {ys, ■ ■ ■ ,yd) and sufficiently small e, the 
set {(^1,2/2): yf + yi < < 0} is path-connected (topologically, it is an annulus). 

Hence C'V(O) U ^^(O) is path-connected. 

Finally, if C3 = and A2C4 < 0, then either the set {yi : V{y) < 0} for fixed (7/2, ■ ■ ■ , yd) 
or the set {y2 ■ V{y) < 0} for fixed (yi, 7/3, ... , y^) is not path-connected, so that the valley 
of is locally split into two disconnected components, joined at the origin. □ 

Remark 2.10. The normal-form transformation x x + g{x) can also be applied when 
Ai 7^ 0. The result is exactly the same normal form as in (j2.13p . except that there is an 
additional term ^AiyJ. This observation is useful as it allows to study the system with a 
unique transformation of variables in a full neighbourhood of the bifurcation point. 

Remark 2.11. One easily checks that if V is of class C, one can construct higher-order 
normal forms by eliminating all terms which are not of the form y\ with k ^ r. In other 
words, there exists a polynomial g{y) with terms of degree between 2 and r — 1 such that 

1 

Viy + giy)) = - ^ X^yf + C^yf + Qy^ + • • • + dyl + o(||y||^) . (2.17) 

i=2 

In general, however, there is no simple expression of the coefficients of the normal form in 
terms of the original Taylor coefficients of V. Note that for z = to be a saddle, the first 
index q such that Cq ^ has to be even with X2Cq < 0. 
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2.4 Singularities of codimension 2 

In this subsection, we shall assume that the potential V is of class C for some r ^ 3 
and that z is a stationary point of V with the Hessian \/'^V{z) having two vanishing 
eigenvalues. We may assume z = and V{z) = 0. According to Proposition 12. 6^ there are 
two cases to be considered: 

1. Ai < 0, A2 = A3 = and < A4 < • • • ^ A^. 

2. Ai = A2 = and < A3 ^ • • • ^ A^. 

For the purpose of this subsection, it will be convenient to relabel the first three eigenvalues 
in such a way that Ai = A2 = 0, while A3 7^ can be positive or negative. We choose 
a basis in which V^F(z) = diag(0, 0, A3, . . . , A^). The potential thus admits a Taylor 
expansion of the form 

1 

V{x) = -'^Xixf + ^ VijkXiXjXk+ ^ VijkiXiXjXkXi + . . . . (2.18) 

i=3 l^i^j^k^d l^i^j^k^l^d 

The theory of normal forms immediately yields the following result. 

Proposition 2.12. There exists a polynomial change of variables x = y + g{y), where g 
is a polynomial with terms of degree between 2 and r, such that 

-, d r 

V{y + g{y)) = -J2 W + E ^^(^1' 2^2) + odlyllD , (2.19) 

i=3 k=3 

where each Vk is a homogeneous polynomial of order k, i.e., 
V3{yi,y2) = Vuiyf + Vii2y?y2 + ^i22yiyi + ^222^2 > 

^4(^1,^2) = Viiiiyf + Viu2yfy2 + ^ii22y?yi + Vi222yiyl + ^2222?/! , (2.20) 

and similarly for the higher- order terms. 

The proof uses standard normal form theory and follows along the lines of the proof 
of Proposition 12.81 given in Appendix |X1 Therefore, we refrain from giving its details 
here. Let us remark that if V is of class C"^^, then (j2.19p holds with 0(11^112) replaced by 
0{\\y\\l^'). 

We can again apply the result to derive an easy to verify necessary condition for a point 
z to be a saddle. Let p be the smallest k such that 14 is not identically zero. Generically, 
we will have p = 3, but other values of p are quite possible, for instance due to symmetries 
(see Section [4.31 for examples). We call discriminant the polynomial! 

Mt) = 1^^^*' = + Vi...utP-' + ■■■ + V2...22 if / , 

\Vp{l,t)=V2...22t'P + V2...2ltP-^ + --- + Vl...ll ifyi...ll=0. 

The following corollary as well as Table [1] provide necessary conditions for z to be a 
saddle, expressed in terms of A(f) and the sign of A3. 



■^The roles of 3/1 and 2/2 being interchangeable, both definitions in (|2.2f |l are equivalent via the trans- 
formation t i-^ 1/t and multiplication by . We choose to make this distinction in order to avoid having 
to introduce "roots at infinity". Equivalently, one could work on the projective line RP^. 
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A3 < 


A3 > 


All roots of A(i) real and simple 
No real root of A(i) and A(i) > 
No real root of A(i) and A(i) < 


Not a saddle 
Saddle 
Not a saddle 


Saddle 

Not a saddle (local minimum) 
Not a saddle 



Table 1. Classification of stationary points with double-zero eigenvalue with the help of 
the discriminant A{t) and the sign of A3. 



Corollary 2.13. 

1. Assume d > 2 and A3 < 0. Then the point z is 

• not a saddle if the discriminant has one or several real roots, all of them simple; 

• a saddle if the discriminant has no real root and is positive; 

• not a saddle if the discriminant has no real root and is negative. 

2. Assume d = 2 or A3 > 0. Then the point z is 

• a saddle if the discriminant has one or several real roots, all of them simple; 

• not a saddle ( a local minimum ) if the discriminant has no real root and is positive; 

• not a saddle if the discriminant has no real root and is negative. 

Proof: In order to determine the open valley OV{z) of z = 0, we first look for solutions 
of V{yi,y2, 0, . . . , 0) = near the origin. Such solutions are necessarily of the form yi = 
y2{t + R{y2)) for some t, where R{y2) goes to zero continuously as y2 — 0. Plugging in, 
dividing by y^ and setting y2 = 0, one sees that t must be a root of the discriminant. 
Applying the implicit function theorem to the pair (^2, R), one finds that there is a unique 
R if this root is simple. Thus whenever the discriminant has real roots, all of them simple, 
V{yi,y2, 0, . . . , 0) changes sign in a neighbourhood of the origin, the regions of constant 
sign being shaped like sectors. If A3 < 0, we have to distinguish between three cases: 

1. The discriminant has simple real roots. In the plane ^3 = • • • = yd = 0, there are 
several disconnected regions in which V is negative. However these regions merge 
when ys becomes nonzero. Hence all regions can be connected by a path leaving the 
plane ys = ■ ■ ■ = y^ = 0, so that the origin is not a saddle. 

2. The discriminant has no real roots and is positive. Then for each fixed {yi,y2,y4, . . .), 
the set of 7/3 such that V is negative is not path-connected, and consequently the set 
{y: V{y) < 0} cannot be path-connected either. By adding the origin, the latter set 
becomes path-connected, showing that is indeed a saddle. 

3. The discriminant has no real roots and is negative. Then either d = 3 and the origin is 
a local maximum, or d > 3 and the set of (7/4, . . . , y^) for which V is negative is path- 
connected for each fixed (yi, y2) 2/3)- Thus any two points in the open valley close to the 
origin can be connected (connect both endpoints to points in the set t/4 = ■ ■ • = = 
by a path with constant (yi, 2/2,^/3); and then connect the two paths within the set 
1/4 = • • • = y^i = 0). Thus the origin cannot be a saddle. 

The proofs are analogous in the case A3 > 0. □ 

We note that if A3 < and the degree p of the discriminant is odd, then the origin is 
usually not a saddle. For the sake of brevity, we do not discuss here cases in which the 
discriminant has nonsimple roots, because then the behaviour depends on higher-order 
terms in the Taylor expansion and there is a large number of cases to distinguish. 
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2.5 Singularities of higher codimension 

Generalisation to nonquadratic saddles of higher codimension is now quite obvious. If the 
potential V is of class C and the first q eigenvalues Ai, . . . , Ag of the Hessian \/'^V{z) are 
equal to zero, the normal form can be written as 

d r 

V{y + giy)) = - A^y? + J] ^^^(yi, • • • + o(||y||^) , (2.22) 

i=q+l fc=3 

where again each 14 is a homogeneous polynomial of degree k. Let p denote the smallest 
k such that is not identically zero. Then the role of the discriminant is now played by 
its analogue 

A{ti, . . . ,tg_i) = Vp{ti,t2,...,tg.i,l) , (2.23) 

and the classification is analogous to the one in Tabled] (with Ag+i instead of A3). Equiv- 
alently, one can study the sign of 1^ on a sphere of constant radius. 



3 First-passage times for nonquadratic saddles 
3.1 Some potential theory 

Let {xt)t be the solution of the stochastic differential equation (jl.ip 
set ^ C M*^, we denote by ta = inf{t > 0: G A} the first-hitting 
A,B CR'^, the quantity 

hA,B{x) =r'^{TA < tb} 
is known to satisfy the boundary value problem 

'LhA,B{x) = forxeiAuBY, 
< hA,Bix) = 1 for X e A , (3.2) 

hA,B{x) =0 for X G B , 

where L = eA — (V1^(-),V) is the infinitesimal generator of the diffusion {xt)t- By 
analogy with the electrical potential created between two conductors at potentials 1 and 
0, respectively, /ia,_b is called the equilibrium potential of A and B. More generally, one 
can define an equilibrium potential h^^, defined by a boundary value problem similar 
to (j3.2p . but with Lh\^ = Xh\^. However, we will not need this generalisation here. 

The capacity of the sets A and B is again defined in analogy with electrostatics as the 
total charge accumulated on one conductor of a capacitor, for unit potential difference. 
The most useful expression for our purpose is the integral, or Dirichlet form, 

cap^(i?) = e [ e-^(^)/^||V/iA,B(x)||i dx =: <I>(^uBr (/ia,b) • (3.3) 
J{AuBy 

We will use the fact that the equilibrium potential hA,B minimises the Dirichlet form 

cap^(B)= inf $(Aui?)-W- (3-4) 

^AU subsets of R"* we consider from now on will be assumed to be regular, that is, their complement is 
a region with continuously differentiable boundary. 



. Given a measurable 
time of A. For set^ 

(3.1) 
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Here 'Ha,b is the space of twice weakly differentiable functions, whose derivatives up to 
order 2 are in L^, and which satisfy the boundary conditions in (|3.2p . 

Proposition 6.1 in |BEGK04] shows (under some assumptions which can be relaxed 
to suit our situation) that if x is a (quadratic) local minimum of the potential, then the 
expected first-hitting time of a set B is given by 

e-^(^)/^/iB.(.),B(2/) dy 
E^{ts} = '-^ jj^. . (3.5) 

The numerator can be estimated by the Laplace method, using some rough a priori es- 
timates on the equilibrium potential hj^^(x),B (which is close to 1 in a neighbourhood of 
X, and negligibly small elsewhere). In the generic situation where x is indeed a quadratic 
local minimum, and the saddle z forms the gate from x to B, it is known that 

/ e-^(^)/^ hB M B{y) dy = —^^ ^ e-^(")/^ [l + 0{e^l^\\oge\)\ , (3.6) 

Jb^ B.W,Bvyy ydit(V2yR) ^ V I s i;j , y ) 

cf. [BEGKOll Equation (6.13)]. Thus, the crucial quantity to be computed is the capacity 
in the denominator. In the simplest case of a quadratic saddle z whose Hessian has 
eigenvalues Ai < < A2 ^ • • • ^ A^, one finds 



cf. [BEGKOil Theorem 5.1], which implies the standard Eyring-Kramers formula (II. 3p . 

In the sequel, we shall thus estimate the capacity in cases where the gate between 
A = Be{x) and i? is a non-quadratic saddle. Roughly speaking, the central result, which 
is proved in Section [5l states that if the normal form of the saddle is of the type 

1 

V{y) = -u^{yi) + U2{y2,...,y,) + - Yl A,y| + 0(||y||^+^) , (3.8) 



2 . 

j=g+l 

where the functions ui and U2 satisfy appropriate growth conditions, then 

g-«2(2/2,...,2/,)/e dy2...dy 



c^PAiB) = e"-^^^^^^ II ,l^[l + R{e)], (3.9) 




for certain Si = 5i{e) > 0, 82 = S2{£) > 0. The remainder R{e) goes to zero as e ^ 0, 
with a speed depending on ui and U2- Once this result is established, the computation of 
capacities is reduced to the computation of the integrals in (13. 9p . 



Remark 3.1. If x is a nonquadratic local minimum, the integral in (j3.6p can also be 
estimated easily by standard Laplace asymptotics. Hence the extension of the Eyring- 
Kramers formula to flat local minima is straightforward, and the real difficulty lies in the 
effect of flat saddles. 
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3.2 Transition times for codimension 1 singular saddles 

We assume in this section that the potential is of class at least, as this allows a better 
control of the error terms. Consider first the case of a saddle z such that the Hessian 
matrix V'^V{z) has eigenvalues Ai = < A2 ^ A3 ^ • • • ^ A^. In this case, the unstable 
direction at the saddle is non-quadratic while all stable directions are quadratic. According 
to Corollary 12.91 in the most generic case the potential admits a normal form 

d 

V{y) = -C,y\ + - \,y] + 0(||y||^) (3.10) 
i=2 

with C4 > 0, i.e., the unstable direction is quartic. (Note that the saddle z is at the origin 
of this coordinate system.) 

We are interested in transition times between sets A and B for which the gate G(A, i?) 
consists only of the saddle z. In other words, we assume that any minimal path 7 G 
7-'(A, -B) admits z as unique point of highest altitude. This does not exclude the existence 
of other stationary points in (DV{z), i.e., the potential seen along the path 7 may have 
several local minima and maxima. 

Theorem 3.2. Assume z is a saddle whose normal form satisfies (jS.lOp . Let A and B 
belong to different path- connected components of OV{z) and assume that G{A,B) = {z}. 
Then 



cap^(i^) = ^^f^^'^/^+^/^e~^(^)/^[l + 0(.V4|iog,|5/4)] , (g.n) 

where T denotes the Euler Gamma function. 

The proof is given in Section [5^ In the case of a quadratic local minimum x, combining 
this result with Estimate (j3.6p immediately yields the following result on first-hitting times. 



Corollary 3.3. Assume z is a saddle whose normal form satisfies (jS.lOp . Let O he one 
of the path- connected components of OV{z), and suppose that the minimum of V in O is 
attained at a unique point x and is quadratic. Let B belong to a different path- connected 
component of OV{z), with G{{x},B) = {z}. Then the expected first-hitting time of B 
satisfies 



Note in particular that unlike in the case of a quadratic saddle, the subexponential 
asymptotics depends on e to leading order, namely proportionally to e~^/^. 

Remark 3.4. 

1. If the gate G{A,B) contains several isolated saddles, the capacity is obtained simply 
by adding the contributions of each individual saddle. In other words, just as in 
electrostatics, for capacitors in parallel the equivalent capacity is obtained by adding 
the capacities of individual capacitors. 
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2. We can extend this result to even flatter unstable directions. Assume that the potential 
V is of class 2p + 1 for some p ^ 2, and that the normal form at the origin reads 



1 

V{y) = -C2,y? + \,yj + , (3.13) 



i=2 



with C2p > 0. Then a completely analogous proof shows that (|3.1ip is to be replaced 
by 



^^P^^^) = ^ 0(.^/^^|loge|(2^+i)/2^)] . (3.14) 

Consequently, if the assumptions of Corollary 13.31 on the minimum x and the set B 
are satisfied, then 



E^Itb] = lilZgP) / 2^A2 . . . A, ^lV(z)-Vi.)]/e ri + 0(gl/2p|log g|(2p+l)/2p)1 _ 

V "i^^l^ '^(^)) 

(3.15) 

Note that the subexponential prefactor of the expected first-hitting time now behaves 
like e-(p-i)/2p. As p varies from 1 to oo, i.e., as the unstable direction becomes flatter 
and flatter, the prefactor 's dependence on e changes from order 1 to order 

Consider next the case of a saddle z such that the Hessian matrix \/'^V{z) has eigenval- 
ues Ai < A2 = < A3 ^ • • • ^ Arf. In this case, all directions, whether stable or unstable, 
are quadratic but for one of the stable directions. According to Corollary 12. 9^ in the most 
generic case the potential admits a normal form 



V{y) = -^|Ai|y? + C,yl + ^ ^ X,y] + 0{\\y\\l) (3.16) 

i=3 

with C4 > 0, i.e., the non-quadratic stable direction is quartic. 

Theorem 3.5. Assume z is a saddle whose normal form satisfies ()3.16p . Let A and B 

belong to different path- connected components of OV{z) and assume that G{A,B) = {z}. 
Then 



Cap^(B) = £(1/^ '[All g^/2-l/4^-n^)/en ^0(gl/4|ipg^|5/4)1 ^^ -^^^ 

2Cy^ V ^3---Ad 

The proof is given in Section [531 In the case of a quadratic local minimum x, combining 
this result with Estimate (13. 6p immediately yields the following result on first-hitting times. 



Corollary 3.6. Assume z is a saddle whose normal form satisfies (j3.16p . Let O be one 
of the path- connected components of OV{z), and suppose that the minimum of V in O is 
attained at a unique point x and is quadratic. Let B belong to a different path- connected 
component of OV{z), with G({x},i?) = {z}. Then the expected first-hitting time of B 
satisfies 



- ^/SISS^"'^ + 0(s'/*|.o,e|V.)] . (3.18) 
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Note again the e-dependence of the prefactor, which is now proportional to e^/^ to 
leading order. A similar result is easily obtained in the case of the leading term in the 
normal form having order y'2' for some 2. In particular, the prefactor of the transition 
time then has leading order e^P~^)l'^P: 

(3.19) 

As p varies from 1 to c«, i.e., as the non-quadratic stable direction becomes flatter and 
flatter, the prefactor's dependence on e changes from order 1 to order ^/£. 

3.3 Transition times for higher-codimension singular saddles 

We assume in this section that the potential is at least of class . Consider the case of a 
saddle z such that the Hessian matrix V^y(z) has eigenvalues 

Ai < = A2 = A3 < A4 ^ • • • ^ Ad . (3.20) 

In this case, the unstable direction is quadratic, while two of the stable directions are non- 
quadratic. Proposition 12.121 shows that near the saddle, the potential admits a normal 
form 

y{v) = -\\Myl + ^3(y2,y3) + ^4(^2,2/3) + \Y. \y] + , (3.2i) 

with V3 and V4 homogeneous polynomials of degree 3 and 4, respectively. If V3 does not 
vanish identically, Corollary 12.131 shows that z is typically not a saddle. We assume thus 
that V3 is identically zero, and, again in view of Corollarv 12.131 that the discriminant 

A(t) = ^2222*^ + ^2223*^ + ^^2233*^ + '^^2333* + ^3333 (3.22) 

has no real roots, and is positive with F2222 > 0. It is convenient to introduce polar 
coordinates, writing 

V4(r cos r sin (/?) = r'^k{Lp) , (3.23) 

where we may assume that k{Lp) is bounded above and below by strictly positive constants 
Kj^ ^ Then we have the following result, which is proved in Section [5.31 

Theorem 3.7. Assume z is a saddle whose normal form satisfies (j3.2ip . with V3 = 
and V4 > 0. Suppose, the discriminant A(t) has no real roots and satisfies V2222 > 0. 
Let A and B belong to different path- connected components of OV{z) and assume that 
G{A,B) = {z}. Then 

(3.24) 

Corollary 3.8. Assume z is a saddle whose normal form satisfies (j3.2ip . Let O he one 
of the path- connected components of OV{z), and assume that the minimum of V in O is 
reached at a unique point x, which is quadratic. Let B belong to a different path- connected 
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component of OV{z), with G{{x},B) = {z}. Then the expected first-hitting time of B 
satisfies 



A:(v3)i/2 

(3.25) 







The prefactor is now proportional to e^^^ instead of being proportional to £^^^, which 
is explained by the presence of two vanishing eigenvalues. 

Remark 3.9. This result admits two straightforward generalisations to less generic situ- 
ations: 

1. If the potential is of class C'^p^^ and the first nonvanishing coefficient of the normal 
form has even degree 2p, p ^ 2, then a completely analogous proof shows that 



2p \pj Jo A;((/p)Vpy A4...Arf 

X [l + 0(eV2P|iog^|{2p+i)/2p)] _ (3^26) 

Consequently, if the assumptions of Corollary 13.81 on the minimum x and the set B 
are satisfied, then 



K^/ro\ - 1 / (27r)4A4 ■ ■ ■ Ad jv(,)_v(^)y, 

\p) Jo HvV^p 

X [l + 0(ei/2p|iQg^|(2p+i)/2p)j _ (3^27) 

2. If all eigenvalues from A2 to Xg are equal to zero, for some g ^ 4 and the first nonvan- 
ishing coefficient of the normal form has even degree 2p, p ^ 2, then the prefactor of 
the capacity has order £°'/2-('?"i){p-i)/2p^ ^j^j involves a {q — 2)-dimensional integral 
over the angular part of the leading term in the normal form. 

The other important codimension-two singularity occurs for a saddle z such that the 
Hessian matrix \/'^V{z) has eigenvalues 

= Ai = A2 < A3 ^ • • • ^ Arf . (3.28) 

In this case. Corollary 12.131 states that z is a saddle when the discriminant of the normal 
form has one or more real roots, all of them simple. As a consequence, there can be more 
than two valleys meeting at the saddle. This actually induces a serious difficulty for the 
estimation of the capacity. The reason is that for this estimation, one needs an a priori 
bound on the equilibrium potential hA_,B in the valleys, some distance away from the 
saddle. In cases with only two valleys, hy^^s is very close to 1 in the valley containing A, 
and very close to in the valley containing B. When there are additional valleys, however, 
one would first have to obtain an a priori estimate on the value of hA_,B in these valleys, 
which is not at all straightforward, except perhaps in situations involving symmetries. We 
will not discuss this case here. 
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(a) (b) (c) 




Figure 3. Saddles and open valleys of the normal- form potential (|4.2p . in a two- 
dimensional case, (a) for A2 > 0, (b) for A2 = and (c) for A2 < 0. The system 
undergoes a transversal pitchfork bifurcation at A2 = 0. 



Arguably, the study of singular saddles satisfying (j3.28p is less important, because 
they are less stable against perturbations of the potential. Namely, the flatness of the 
potential around the saddle implies that like for the saddle-node bifurcation, there exist 
perturbations of the potential that do no longer admit a saddle close by. As a consequence, 
there is no potential barrier creating metastability for these perturbations. 

4 Bifurcations 

While the results in the previous section describe the situation for nonquadratic saddles, 
that is, at a bifurcation point, they do not incorporate the transition from quadratic to 
nonquadratic saddles. In order to complete the picture, we now give a description of 
the metastable behaviour in a full neighbourhood of a bifurcation point of a parameter- 
dependent potential. We will always assume that V is of class C^. 

We shall discuss a few typical examples of bifurcations, which we will illustrate on the 
model potential 

N N 

y^(x) = ^C/(x,) + J^(x,-x,+i)2, (4.1) 

i=l 1=1 

introduced in [BFGOTaj . Here Xj denotes the position of a particle attached to site i of the 
lattice Z /NZ, , U{xi) = ^xf — ^xf is a local double-well potential acting on that particle, 
and the second sum describes a harmonic ferromagnetic interaction between neighbouring 
particles (with the identification x^r+i = 2:1). Indeed, for N = 2 and 7 = 1/2, the origin 
is a nonquadratic saddle of codimension 1 of the potential (j4.ip . while for all ^ 3, the 
origin is a nonquadratic saddle of codimension 2 when 7 = (2 sin^(7r/A^))^^. 

4.1 Transversal symmetric pitchfork bifurcation 

Let us assume that the potential V depends continuously on a parameter 7, and that for 
7 = 7*, 2; = is a nonquadratic saddle of V, with normal form (j3.16p . A symmetric 
pitchfork bifurcation occurs when for 7 near 7*, the normal form has the expression 

Viy) = ^Ai(7)y? + ^A2(7)yi + ^^4(7)^2' + ^ E ^^(^^y' + ^(H^H^) ' (4-2) 

where A2(7*) = 0, while Ai(7*) < 0, (74(7*) > 0, and similarly for the other quantities. We 
assume here that V is even in 1/2, which is the most common situation in which pitchfork 
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Figure 4. Level lines of the potential (|46)) . (a) for 7 > 1/2, (b) for 1/2 > 7 > 1/3 
and (c) for 1/3 > 7 > 0. A transversal pitchfork bifurcation occurs at 7 = 1/2, and a 
longitudinal pitchfork bifurcation occurs at 7 — 1/3. 



bifurcations are observed. For simplicity, we shall usually refrain from indicating the 7- 
dependence of the eigenvalues in the sequel. All quantities except A2 are assumed to be 
bounded away from zero as 7 varies. 

When A2 > 0, 2: = is a quadratic saddle. When A2 < 0, z = is no longer a saddle 
(the origin then having a two-dimensional unstable manifold), but there exist two saddles 
z± with coordinates 

z± = (0, ±^1^21/4^4 + 0(A2), 0, . . . , 0) + 0(Al) (4.3) 

(Figure [3|). Let us denote the eigenvalues of V'^V{z±) by /ii, . . . , Hd- In fact, for A2 < 
we have 

^2 = -2A2 + 0(|A2|'/2) , 

/i, = A,+0(|A2|3/') for jG {1,3,..., 4. (4.4) 

Finally, the value of the potential on the saddles z± satisfies 

V{z+) = V{z.) = Viz) - -Jr + 0{\X2\'/') . (4.5) 



Example 4.1. For N = 2 particles, the potential (j4.ip reads 

V{X^,X2) = U{xi) + U{X2) + ^{Xi - X2? . (4.6) 

Performing a rotation by 7r/4 yields the equivalent potential 

v{yi.y2) = -\yl- ^-^yl + \{yt + ^yhl + y\) , (4.7) 

which immediately shows that the origin (0,0) is a stationary point with Ai(7) = — 1 and 
A2(7) = — (1 — 27). For 7 > 7* = 1/2, the origin is thus a quadratic saddle, at "altitude" 0. 
It serves as a gate between the local minima located at y = (ibl/\/2, 0). As 7 decreases 
below 7*, two new saddles appear at y = (0, ±y^2(1 - 27) (cf. Fi gured] which shows the 
potential's level lines in the original variables (xi,X2)). They have a positive eigenvalue 
/^2(7) = 2(27 — 1), and the "altitude" —^(1 — 27)^. There is thus a pitchfork bifurcation 
at 7 = 1/2. Note that another pitchfork bifurcation, affecting the new saddles, occurs at 
7 = 1/3. 
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Figure 5. The functions 'i'±{ct), shown on a hnear and on a logarithmic scale. 



Our main result is the following sharp estimate of the capacity. 



Theorem 4.2. Assume the normal form at z = satisfies ()4.2p . Let A and B belong to 
different path- connected components ofOV{z) (respectively ofOV{z±) if X2 <0)- Assume 
further that G{A, B) = {z] (resp. G{A, B) = {z-,z+} if A2 < Oj. Then for A2 > 0, 



cap^(s) = y 

while for A2 < 0, 
cap^(B) = 



(27r)rf-2|Ai| 



[A2 + (2eC4)V2]A3...A, ^V(2eC4)i/2 



A2 



e'^/'e-^(^)/ni + i?+(e,A2)] , 



(4.8) 



(2vr)'^-2|^i| 



V [/.2 + (2eC74)V>3...A^d l(2eC4)V2f [i + ^-(e,^2jJ . 

(4.9) 

The functions and are hounded above and below uniformly on M + . They admit 
the explicit expressions 



a(l + a 



7ra(l + a) _q,2 



/64 



32 



(a' 



(4.10) 



where i^i/4 and I±i/4 denote the modified Bessel functions of the second and first kind, 
respectively. In particular, 



and 



lim ^j^{a) = l, lim ^^-(a) = 2 



r(i/4) 

lim ^'+(a) = lim ^'_(a) = ' ~ 0.8600 . 

a^o Q^o 2^' v/vr 



Finally, the error terms satisfy 

\R±ie,X)\^C 



e|loge| 



max{|A|,(e|loge|)V2} 



1/2 



(4.11) 
(4.12) 

(4.13) 
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Figure 6. The prefactor of the expected transition time near a pitchfork bifurcation, as 
a function of the bifurcation parameter A2, shown for three different values of e. (To be 
precise, we show the function A2 ^ ■\/A2+^^/^'+(A2/£^^^) for A2 > and the function 
A2 ^ V-2A2 + ei/V*-(-2A2/ei/2) for A2 < 0.) 



The functions \I'±(a) are universal, in the sense that they will be the same for all 
symmetric pitchfork bifurcations, regardless of the details of the system. They are shown 
in Figure[5j Note in particular that they are not monotonous, but both admit a maximum. 



Corollary 4.3. Assume the normal form at z = satisfies (j4.2p . 

• If ^2> 0, assume that the minimum of V in one of the path- connected components of 
OV{z) is reached at a unique point x, which is quadratic. Let B belong to a different 
path- connected component of OV{z), with G{{x},B) = {z}. Then the expected first- 
hitting time of B satisfies 



[A2 + (2eC4)i/2]A3 



A. 



,[V(z)~V(x)]/e 



|Ai|det(V2y(x)) ^+(A2/(2eC4)V2) 



[l + R+{e,\2)\ . (4.14) 



If ^2 < 0, assume that the minimum of V in one of the path- connected components of 
OV{z-^-) = OV{z-) is reached at a unique point x, which is quadratic. Let B belong to 
a different path- connected component of OV{z±), with G{{x},B) = {z^,Z-}. Then 
the expected first-hitting time of B satisfies 



[lJ-2 



(2£C4)V2]^3, 



^J'd 



^[V{z±)-V{x)]/s 



|A^i|det(W(x)) ^'_(/i2/(2eC4)i/2) 



[l + i?_(e,^2)] . (4.15) 



When A2 is bounded away from zero, the expression (j4.14p reduces to the usual Eyring- 
Kramers formula (jl.3p . When A2 — > 0, it converges to the limiting expression (j3.18p . The 
function ^_(_ controls the crossover between the two regimes, which takes place when A2 
is of order e^/^. In fact, when A2 <C e^^'^, there is a saturation effect, in the sense that 
the system behaves as if the curvature of the potential were bounded below by (2eC4)^/^. 
Similar remarks apply to the expression (|4.15p . the only difference being a factor 1/2 in 
the prefactor when /X2 is bounded away from (cf. ()4.1ip l. which is due to the fact that 
the gate between x and B then contains two saddles. 

The A2-dependence of the prefactor is shown in Figure [H It results from the combined 
effect of the term under the square root and the factors ^'-i-. Note in particular that the 
minimal value of the prefactor is located at a negative value of A2, which can be shown to 
be of order e^/^. 
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(a) (b) (c) 




Figure 7. Saddles and open valleys of the normal-form potential (|4.16p . (a) for Ai < 0, 
(b) for Ai = and (c) for Ai > 0. The system undergoes a longitudinal pitchfork 
bifurcation at Ai = 0. 



4.2 Longitudinal symmetric pitchfork bifurcation 

Consider now the case where for 7 = 7*, z = is a nonquadratic saddle of V, with normal 
form (|3.10p . Then a slightly different variant of symmetric pitchfork bifurcation occurs, 
when for 7 near 7*, the normal form has the expression 

V{y) = lxi{i)yl - C,{^)y\ + \Y1 W?^' + ^dl^^H^) ' (4-16) 

i=2 

where Ai(7*) = 0, while (74(7*) > and Aj(7*) > for all j ^ 2. We assume here 
that V is even in yi, which is the most common situation in which pitchfork bifurcations 
are observed. As before, we no longer indicate the 7-dependence of eigenvalues, and all 
quantities except Ai are assumed to be bounded away from zero. 

When Ai < 0, z = is a quadratic saddle. When Ai > 0, 2: = is a local minimum, 
but there exist two saddles z± with coordinates 

z± = (± v^AlTi^ + 0(Ai), 0, . . . , 0) + 0{\l) . (4.17) 

The open valleys of and z_ share a path-connected component which we denote by Oq, 
while we denote their other components by 0±, cf. Figure [71 Let us denote the eigenvalues 
of V^y(z-i-) by /Ui, . . . , /xrf. For Ai > we have 

//I = -2Ai + 0(Af , 

lij=\j + 0{\^^^) for j E {2,... ,4 . (4.18) 

Finally, the value of the potential on the saddles z± satisfies 

V{z+) = V{z.) = V{z) + + 0{\T) . (4.19) 

Such a bifurcation occurs, for instance, in Example 14.11 for 7 = 1/3. Then both 
saddles z^ and Zj^ have to be crossed on any minimal path between the global minima, as 
in Figure [1};. 

We can now state a sharp estimate of the capacity in this situation, which is proved 
in Section 15.41 

Theorem 4.4. Assume the normal form at z = satisfies (14.160 . Let A and B belong to 
the pathconnected components O- and 0+ respectively. Assume further that G{A, B) = 
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Figure 8. The prefactor of the expected transition time near a longitudinal pitchfork 
bifurcation, as a function of the bifurcation parameter Ai, shown for three different values 
of £. (To be precise, we show the function Ai ^ ^^{\\r\le^^'^)l ^\\x\^ e^l"^ for Ai < 
and the function Ai *_(2Ai/ei/2)/^2Ai + for Ai > 0.) 



{z} (resp. that {z-} and {z_|_} both form a gate between A and B if Xi > 0). Then for 
Ai <0, 



cap^(S) 



/ ATTTA^ ^+(|Ai|/(2.C4)V2) ' + ' 

(4.20) 

while for Ai > 0, 



(27r)^-2[|/ii| + (2£C4)V^] g-^/^ -y(.+)/er-, , o r , m 

^^^^^'') = V ^-(|..|/(2.C4)V^) ^ + |,,|)]. 

(4.21) 

T/ie functions anc? anrf i/ie remainders R± are the same as in (j4.10p and (|4.13p . 

When Ai ^ e^/^, the function in (j4.2ip is close to 2, so that the total capacity 
equals half the capacity associated with each saddle Z- and z^. As in electrostatics, when 
two capacitors are set up in series, the inverse of their equivalent capacity is thus equal to 
the sum of the inverses of the individual capacities. 

Corollary 4.5. Assume the normal form at z = satisfies (14.161) . 

• If Xi < 0, assume that the minimum of V in one of the path- connected components of 
OV{z) is reached at a unique point x, which is quadratic. Let B belong to a different 
path- connected component of OV{z), with G({x},i3) = {z}. Then the expected first- 
hitting time of B satisfies 



^2---Xd ( I All \ \y(z)-V{x)'\le 



W'{tb) 27tJ [i^^i + (2eC74)V2]det(V2y(x))'^+V(2eC4)i/2 ) ^ 



[l + i?+(e,|Ai|)] . (4.22) 



• If Xi > 0, assume that the minimum of V in the path- connected component O of 
OV{z-) is reached at a unique point x, which is quadratic. Let B belong to the path- 
connected component 0+ of OV{zj^), with G{{x},B) = {zj^} or {z-}. Assume finally 



24 



that V{x) is strictly smaller than the minimum of V in the common component Oq of 
OV{z-) and OV{zj^). Then the expected first-hitting time of B satisfies 



E^TrI = 2-K I 1^2- --f^d ~ / l^ll ^ ^[V{z±)~V{x)]/e 

^ ^ V [lml + (2e<^4)i/2]det(V2F(x)) ~\{2eC^y/y 

X [l + i?_(e,|^i|)] . (4.23) 

For Ai = 0, this result reduces to Corollary 13.31 For negative Ai of order 1, one 
recovers the usual Eyring-Kramers formula, while for positive Ai of order 1, one obtains 
the Eyring-Kramers formula with an extra factor 2, which is due to the fact that two 
saddles have to be crossed in a row. Note that when applying this result to Example 14.11 
near 7 = 1/3, all expected transition times have to be divided by 2, because there are two 
equivalent paths from one global minimum to the other one. 



4.3 Bifurcations with double-zero eigenvalue 

We consider finally the case where the 7-dependent potential V admits, for 7 = 7*, z = 
as a nonquadratic saddle with two vanishing eigenvalues. Then the normal form is given 
by ()3.2ip with V3 = 0. This singularity having codimension 2, there are several different 
ways to perturb it, and actually two parameters are needed to describe all of them. We 
restrict our attention to the particular perturbation (which is generic, e.g., in cases with 
i^iv-symmetry. A*" ^ 3) 

V{y) = ^Ai(7)2/? + ^A2(7)(yi + yl) + V4{y2,y3-n) + ^ E ^^^2/1 + ^dl^Hs) ' (4-24) 

i=4 

where V4(y2,y3;7) is of the form (j2.20p with positive discriminant, satisfying Vim > 0. 
Here A2(7*) = 0, while Ai(7*) < 0, A4(7*) > 0, and so on. Again, we no longer indicate 
the 7-dependence of the eigenvalues in the sequel. All quantities except A2 are assumed 
to be bounded away from zero. 

Example 4.6. Consider the potential (j4.ip for ^ 3. In Fourier variables 

1 ^ 

^'^ = 7^ E ^""''"^"^ ' - L^/2J < ^ ^ LA^/2J , (4.25) 
the potential takes the form 

= 2 XI ^'^^'^^-'^ + 4/V E Zk^Zk^Zk^Zki , (4.26) 

k fci+A:2+fc3+fc4=0 (mod N) 

where rjk = — 1 + 27 sin^(A;7r/A^). For 7 > (2 sin'^(7r/A^))^-'^, all r]^ except r]o = —1 are 
positive, and thus the origin is a quadratic saddle. It forms the gate between the two 
global minima of the potential, given (in original variables) by = ±(1, . . . , 1). As 7 
approaches (2 sin^(7r/A'^))~^ from above, the two eigenvalues rj±i go to zero, and the origin 
becomes a singular saddle of codimension 2. In [BFGOTa] it is shown that as 7 decreases 
further, a certain number (which depends on A^) of quadratic saddles emerges from the 
origin, all of the same potential height. The set of all these saddles then forms the gate 
between I~ and . 
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The normal form of (|4.26|) is equivalent to (j4.24p with d = A^, up to a linear transfor- 
mation (one has to work with the real and imaginary part of each Zk instead of the 
pairs {zk,z^k)), and to a relabelling of the coordinates, setting Ai = ?7o = —1 and 

We start by considering saddles with normal form (14.24p for A2 ^ 0, when the origin 
is a saddle, and for slightly negative A2. Indeed, in these cases one can derive a general 
result which does not depend on the details of the nonlinear terms. We define k{ip), as 
in (|3.23p . by V4{r cos ip,r simp; 'y) = r'^k{ip;j). 

Theorem 4.7. Assume the normal form at z = satisfies (|4.24p . Let A and B belong to 
different path- connected components ofOV{z) (respectively of the newly created saddles if 
A2 < Oj. Assume further that the gate G{A,B) between A and B is formed by the saddle 
z in case A2 ^ 0, and by the newly created saddles, otherwise. Then for A2 ^ 0, 



27r 



X e 



''/'e-^(^)/^[l + i?+(£,A2)] , (4.27) 



while for — (e|log ej)-*^/^ ^ A2 < 



X ^-v{z)/e ^ ^_ ^ (4^28) 

The functions 0+ and 0_ are bounded above and below uniformly on M + . They are given 
by 

e+(a) = + a) e"V8 ^(-fj , (4.29) 

Q_{a) = .[l4^\ , (4.30) 



2 V2, 

where <I>(x) = (27r)^-'^/2 e"^^/^ dy denotes the distribution function of the standard 
normal law. In particular, 

hrn e+(a) = l, lim e-(a) = ./f, (4.31) 



and 



lim e+(Q) = lim Q_{a) = J- ~ 0.6267 . (4.32) 



Finally, the error terms satisfy 

e|loge| 



\R±{e,\)\ ^ C 



{|A|,(e|loge|)V2} 



max 



1/2 

(4.33) 



The proof is given in Section (j5.4p . The functions Q±{a), which are shown in Figure El 
are again universal in the sense that they will be the same for all bifurcations admitting 
the normal form (j4.24p . 



26 




20 0.1 1 10 100 

Figure 9. The functions Q±{ct), shown on a hnear and on a logarithmic scale. 



Corollary 4.8. Assume the normal form at z = satisfies (|4.24p . 

• If ^2 > 0, assume that the minimum of V in one of the path- connected components of 
OV{z) is reached at a unique point x, which is quadratic. Let B belong to a different 
path- connected component of OV{z), with = {z}. Then the expected first- 

hitting time of B satisfies 



A4 



A. 



,[V{z)-V{x)ye 



Ai|det(V2y(x)) 1 /-s- e+(A2/(2eA:((^))V2) 



■[l + i?+(e,A2)] . 



27r io As + (2eA:(v9))V2 



dtp 



(4.34) 

The above situation extends to slightly negative X2, that is, — Y^e|loge| ^ A2 ^ 0, 
where 



Xa . . . X„ 



^[Viz)-Vix)]/e 



|Ai|det(V2T/(x)) 1 
2^ 



2tt 



e_(A2/(2eA;(vp))V2 



(2eA:(99))i/2 



dip 



[l + R_ie,X2)] . (4.35) 



For positive A2 of order 1, one recovers the usual Eyring-Kramers law (recall that 
A2 = A3), while for A2 = 0, one recovers Corollarv 13.81 For negative A2, if the function 
k{ip) is nonconstant, the integral over ip in (j4.35p can be evaluated by the Laplace method. 
Extrapolating the result to A2 of order —1, the integral would yield an extra e^^'^ which 
cancels with the the integral's denominator. 



Example 4.9. Consider the potential ()4.26p for = 3 or ^ 5. The resonant terms 
near the bifurcation occurring for Ai = are the monomials proportional to ziZ-i and to 
zfz'^i- In order to obtain the standard normal form (I4.24p . we set z±i = (2/2 i iy^)/^/2 
(the factor \/2 guarantees that the change of variables is isometric) . This yields 



3 



(4.36) 



and thus a constant k{p) = 3/8N. The integrals in (j4.34p and (j4.35p can thus easily be 
computed. For instance, for odd N and A2 = — 1 + 27 sin2(7r/A^) > 0, the prefactor of the 



27 



transition time from I to a neighbourhood of is given by 

\^K...\n-i A2 + (3e/4iV)V^ . ^j, ^^^ (a .7^ 

Vdet V^y(/-) B,(A./(3./4iV) V^) + ^^)] ' ^''''^ 

which admits a limit as A2 ^ 0+. Note that the eigenvalues of V^y(/~) are of the form 
vi = 2, = V2k+i = 2 + 27sin^(A;7r/A^), so that all quantities in (j4.37p are known. 

Example 4.10. Consider now the potential (j4.26p for = 4. Then there are two addi- 
tional resonant terms, namely the monomials proportional to z\ and zl^. Proceding as in 
the previous example, this yields 

^4(2/2,^3) = ^(2/1 + ^1) , (4.38) 

and thus k{Lp) = (3 + cos(4(/?))/32. It is, however, more convenient to keep rectangular 
coordinates instead of using polar coordinates. Then the coordinates y2 and 1/3 separate, 
each one undergoing independently a pitchfork bifurcation. 



The behaviour for negative A2 < — y^e|loge| is determined by the shape of the potential 
in the (2/2, y3)-plaiie. Near the bifurcation point, the potential is sombrero-shaped. As A2 
decreases, however, the sombrero may develop "dips in the rim", and these dips will 
determine the value of the integral over 99. We conclude the discussion by an in-depth 
study of this phenomenon in the case of our model potential. 

Example 4.11. Consider again the potential (j4.26p for = 3 or ^ 5. As seen in 
Example 14.91 above, the quartic term then is rotation-invariant. This, however, is not 
sufficient to determine the behaviour for A2 < — y^e|loge|. 

In fact, we know from symmetry arguments [BFGOTa] that the normal form of the 
potential around the origin has the form 

M N 

y{y) = -^1 Ally? + '-\2r' + E C2,r^' + D2Mr'"' cos(2M^) + ^ A,y| + 0{\\y\\f"+^) , 

q=2 j=4 

(4.39) 

with M = if is odd, and M = N/2 if A^ is even. Here (r, denote polar coordinates 
for (2/2,2/3)- To be precise, the approach presented in [BFGOTa] is based on a centre- 
manifold analysis, but the expression of the potential on the centre manifold is equal to 
the resonant part of the normal form (this is because after the normal-form transformation, 
the centre manifold is ||y||2*^'''^-close to the (7/2, y3)-plane). 

For A2 < 0, V{y) has nontrivial stationary points of the form z* = (0, r* , if* , 0, . . . , 0) -|- 
C(A|), where 

= .lzh + o{\\2\'"^) , sin(2M<^*) = . (4.40) 
V 4C4 

The eigenvalues of the Hessian around these points are of the form 
W = Ai + 0(|A2|) , 

^2 = -2^ir\^^) = ±{2MfD2M[j^) +0{\X2n , 

^fc = Afc + 0(|A2|) , fc = 4,...,Ar. (4.41) 
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Figure 10. A2-dependence of the prefactor of the expected transition time for the model 
system gl]) for iV = 3. 



The points with 112 > are saddles. They correspond to the "dips on the rim of the 
sombrero". In the sequel we will therefore assume fi2 > 0. 

Defining the sets A and B as usual, an analogous computation to the previous ones 
(see Section 15.41 for details) shows that 



cap^(.B) = 2M4 



(27r)'^-2|/.i 



(- 



X £ 



d/2 ^-V{z*)/e 



[l + i?_(e,M2)] , (4.42) 



with i?_(e,/i2) as in ()4.33p . Here 0_ is the function defined in ()4.30p . and 



X(a) =2^/rT^e-"/ol 



a] 



where Iq is the modified Bessel function of the first kind. It satisfies 



lim x(a) = 2 , lim x(a) 



It follows that if i? is a small neighbourhood of I"*", then for A2 < 0, 

' ^[V{z*)-V{x)\/e 



^^^J 2M V 



|/xi|det(V2T/(/-)) 



e. 



^3 



V(8eC4)i/V^V(2M)28eC4, 

X [l + i?_(e,^2)] . (4.45) 



/"2At3 



(4.43) 



(4.44) 



The normal- form coefficient C4 is in fact equal to 3/8N (compare (|4.36p ). which allows to 
check the continuity of the prefactor at A2 = 0. One can identify three different regimes, 
depending on the value of ^3 > (or A2 < 0): 

• For ^ /is ^ ^/e, the arguments of the functions 0- and x are negligible, so that 
the denominator can be approximated by 0_(O)x(O) = Y^7r/2. Furthermore, the 
term (2M)^8eC4 under the square root dominates /X2/U3. This results in a prefactor 
of order e^^'^, compatible with (|4.35p and (j4.37p . In this regime, the potential in the 
(2/2, 2/3)-plaiie is almost flat. 
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• For yje <C ^3 ^ e^/^^, the argument of G_ is very large, so that 0_ can be ap- 
proximated by y^7r/2, while the argument of x is still negligible. The product of 
0_ and X is thus close to \/27r. In this regime, the potential in the (?/2, y3)-plane is 
sombrero-shaped, and very close to rotation-invariant. 

• Finally, for //3 e^/^^, the arguments of both 0_ and x are very large, so that their 
product is close to 1. Furthermore, /X2//3 under the square root dominates (2M)^8eC4. 
In this regime, the potential in the (^2; y3)-pla'ne is sombrero-shaped, but with 2M 
dips of noticeable depth. One recovers the standard Eyring-Kramers formula, with an 
extra prefactor of 1/2M accounting for the fact that the gate consists of 2M saddles 
at equal height. 

The A2-dependence of the prefactor is shown in Figure [10] in the case = 3. In 
general, the prefactor behaves like A2 for A2 ^ \/e and like |A2|^''''^ for A2 ^ —e^^^'^. Note 
that as N increases, the second crossover at e^/*^ occurs later and later, as the system 
becomes closer and closer to being rotation-invariant. 
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5 Proofs 



5.1 Upper bound on the capacity 

We assume that the potential V is of class for some r > 2, and that the origin is 
a stationary point with eigenvalues satisfying Ai ^ < A^+i ^ • • • ^ A(i for some q ^ 2. 
The eigenvalues A2, . . . , Ag lie between Ai and Ag+i, but may have an arbitrary sign. We 
assume that Xq+i ^ • • • ^ A^ are of order 1. In the vicinity of the origin, V admits a 
normal form 

1 

V{y) = -ui{yi)+U2iy2,...,yg) + ^ J2 ^jv] + 0{\\y\\l-^') , (5.1) 

where we may assume that ui and U2 are polynomials of degree less or equal r, and 
increasing at infinity. Suppose that 

• either the origin z* = itself is a saddle, 

• or that there exists a finite number of saddles z* in the vicinity of z = 0. 

Let A and B belong to two different path-connected components of OV{z*) or OV{z*), re- 
spectively; and assume for simplicity that the gate G{A, B) consists of the above-mentioned 
saddles only. The assumptions below will guarantee that the saddles z* arc close to 2; = 0. 
However, the results can easily be extended to more complicated situations with several 
gates which are not necessarily close to the origin, simply by summing the contributions 
of all gates. 

Proposition 5.1 (Upper bound). Assume there exist strictly positive numbers 5i = 5i{e), 
<^2 = ^2 (s) and c'^Q ( independent of e) such that 

""1(2/1) ^ del log e| whenever \yi\ ^ 5i , 

^2(2/2, • • • ,yg) ^ -c<ie|loge| whenever ||(y2, • • • ,yg)||2 ^ (^2 , 

■"2(2/2, 2/q) ^ 2d£|loge| whenever \\{y2,..., yq)\\2 ^2 , (5.2) 



and such that 
Then 



[5i{e) + 52{e)Y^^ = o{e\\oge\) . (5.3) 



j c^«2{j/2,...,J/,)/£ (iy^ . . . dyg ^ . 

cap^(^)^6 "'-^'^ n \r-^[l + Ri{e)\+R2{e), (5.4) 

/ Q-ui{yi)/e j=q+l V J 

J-Si 

where 13$^ (0) = {{y2, ■ ■ ■ , yq) : 2/2 + ' ' ' + 2/g < '^2} ^^'^ error terms satisfy 



Ri{e) ^ 



'c[eV2|iog£|3/2 + e-\5l+^ + <5^+i) + (5i + 82] for r = 2, 
C[£V2|iog£|i/2 + £-1(^5;+! + + 5, + ,52] for r > 2, 

R2{e) ^ Ce^^l^+^Sf ^ Ce^'^+^'H^ (5.5) 



for some constant C > 0. 



31 



Proof: The proof is adapted from the proof of [BEGK041 Theorem 5.1]. 

Recah that the capacity can be computed as the minimal value of the Dirichlet form 
^{AuBy^ cf. (j3.4p . which involves an integration over x. In the vicinity of the saddle, we 
carry out the normal-form transformation of Proposition 12.81 or 12.121 in the integral. This 
can always be done locally, setting x = y + p o g{y), where p is a smooth cut-off function 
which is the identity in a small ball of radius A, and identically zero outside a larger ball 
of radius 2A. Inside the smaller ball, we may thus assume that the potential is given 
by dSI]). 



Let 6 = 2y/ (1 + c)de|loge|. We introduce a set 



d 

C, = [-5i,5i\ X BsM X W [-5j, 5j] , (5.6) 

j=q+l 



where 5i and 82 satisfy (15. 2p and (15. Sp and we choose 

5 



5j = —= , j = q+l,...,d . (5.7) 



^3 

By assumption, making e small enough we can construct a layer of width 26i, separating 
the connected components of the open valley C'V(O), and such that V{y) is strictly positive 
for all y £ iS^yCe. Let Z)_ and denote the connected components of R'^\5£ containing 
A and B respectively. Note that we can choose a radius A for the ball -Ba(O) in which 
we carry out the normal form transformation independently of £ in such a way that V is 
bounded away from zero on Ss H -Ba(O)'^. 

The variational principle (j3.4p implies that it is sufficient to construct a function G 
'Ha,b such that ^(AuB^ih'^) satisfies the upper bound. We choose 



h+iy) 



1 for y G D_ , 

for y G Z?+ , (5.8) 

[f{yi) foiyeCe 



while h'^iy) is arbitrary for y £ S^X Cg, except that we require ||V/i^||2 ^ const /5i. The 
function f{yi) is chosen as the solution of the one-dimensional differential equation 

£r(yi)-|^(yi,o,...,o)/'(yi) = o (5.9) 

with boundary conditions 1 in —61 and in that is, 

'"''e^(*'0'-'0)/^dt 
Si ' 



fiyi) = • (5-10) 



^Vit,0,...,0)/e 

Si 

Inserting into the expression (j3.3p of the capacity, we obtain two non- vanishing terms, 
namely the integrals over 5e\C£ and over C^. The first of these can be bounded as follows. 
For sufficiently small A, Assumptions (|5.2p and (j5.3p imply that for all y £ {Ss\Cs)r\B/\{0), 

^ ^ -d\loge\ - cd\loge\ + 2(1 + c)fi|loge| + 0{£-^[5 + di + 62]"^^) 

>ld\loge\ (5.11) 
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for sufficiently small e. On (5^ \ C^) n B^{Oy, V is bounded away from zero, and by the 
assumption of exponentially tight level sets, the contribution from this part of the integral 
is negligible. It follows that 

, f ^-Viy)/e 2^dy = ©(e^'^/^+iJr') =: Me) • (5-12) 

Js,\c, 0^ 

When calculating the second term which is given by 



^-V{y)/e ^2V{y^,0,...,0)/e 

^cAh^) = e [ e-^(^)/^|/'(yi)|2 dy = e^^^ , (5.13) 

■^^^ e^(^i'O'-'O)/^ dyij 

it is sufficient to carry out the normal form transformation in a smaller ball of radius 2A 
with A = y/dmax{S, 61,62} as -B^(O) contains C^. The Jacobian of the transformation 
thus yields a multiplicative error term 1 + 0{A). By ()5.ip . we have for y G 



Viy)-2Viyl,0,...,0) = u^{y^)+U2{y2,...,yg) + ^ ^ \jy]+0 {[6+61+62]'+^) . (5.14) 

j=q+l 

Hence the numerator in (j5.13p is given by 



e-^i^S'i)/^ dyi [ ^-My2,-,y,)/s . . . dy, f[ e'^^y^^^' dyj 

'[6 + 61+ 62]"-+^ 



j=q+l 
1 + 



(5.15) 



Substituting in (I5.13P we get 



f ^-U2iy2,...,y,)/e ^y^ ^ ^^ ^y^ ^ 

Jb 



^cAh'-)=e'-^^^^ n / e-^^^./^^dy, 



(5A6) 

Using the fact that the Gaussian integrals over yj, j = q + 1, . . . , d, are bounded above by 
Y^27re/Aj, the desired bound (j5.4p follows. □ 



Remark 5.2. Using the conditions (15.21) in order to bound the integrals over yi and 
(2/2, • • • ) Uq), one obtains as a rough a priori bound 

cap^(B) ^ consi^ei-'?/2-('=+i/2)rf[i + +^2(e) . (5.17) 

61 

In applications we will of course obtain much sharper bounds by using explicit expressions 
for lii(yi) and ^2(2/2, • • • 1 yg)i but the above rough bound will be sufficient to obtain a 
lower bound on the capacity, valid without further knowledge of the functions ui and U2- 
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5.2 Lower bound on the capacity 



Before we proceed to deriving a lower bound on the capacity, we need a crude bound 
on the equilibrium potential We obtain such a bound by adapting similar results 

from |BEGK04i Section 4] to the present situation. 

Lemma 5.3. Let A and B he disjoint sets, and let x G (A U BY be such that the ball 
does not intersect AU B. Then there exists a constant C such that 



,V({x},B)/e 



hA,B{x) ^ Ce "'capg^(^)(A)e^ 
Proof: jBEGK04l Proposition 4.3] provides the upper bound 



hAM^) ^ C 



(5.18) 



(5.19) 



so that it suffices to obtain a lower bound for the denominator. This is done as in [BEGK041 
Proposition 4.7] with p = £, cf. in particular Equation (4.26) in that work, which provides 
a lower bound for the capacity in terms of an integral of e^/*" over a critical path from 
X to B. Evaluating the integral by the Laplace method, one gets e^«^'>'-^)/^ as leading 
term, with a multiplicative correction. The only difference is that while Bovier et al 
assume quadratic saddles, which yields a correction of order ^/£, here we do not assume 
anything on the saddles, so that in the worst case the prefactor is constant. This yields 
the bound KW\ . □ 



The capacity capg^(2,)(yl) behaves roug hly like e-V({x\,A)/e ^ ^^^^^ ^j^g bound ([538]) 
is useful whenever V{{x},A) » V{{x},B). This is the case, in particular, when A and 
B belong to different path-connected components of the open valley of a saddle z, and x 
belongs to the same component as B. If, by contrast, x belongs to the same component as 
A, the symmetry hA,B{x) = 1 — hs^Aix) yields a lower bound for the equilibrium potential 
which is close to 1. 

We now consider the same situation as in Section [5Tl Let Si{e), 62(e) and c be the 
constants introduced in Proposition 15.11 

Proposition 5.4 (Lower bound). Let K > max{2d(2 + c), d — g}, and assume that 5i{e) ^ 
£0 Furthermore, suppose there exist strictly positive numbers 61 = (5i(e) and 82 = 62{s) 
such that 



ui{±5i) ^Ke\\oge\ , 
U2{y2,---,yq) ^ Ke\loge\ , 
and such that 

[5i{e)+52{e)Y+^ = o{e\loge\) 

Then 



whenever \\{y2, . . .,yq)\\2 ^ 62 , 



(5.20) 



capj^{B) ^ e 



I e-"2(y2.-.f9)/e d?/2...dy, 



9 d 



Si 



<5i 



-ui{yi)/e 



n 

j=q+l 




[l-Rsie)], (5.21) 



*It actually suffices to have 5i ^ ^k:/2+i/4 d/s ^jjjgjj jg ^ very weak condition satisfied whenever 5i ^ 
for a K > 1 depending on c, d and q. 
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where the remainder i?3(e) satisfies 

\c[ei/2|loge|V2+e-i(^^'+i + ^^'+i) + ^, +<52 + V^] forr>2, 
for some constant C > 0. 



(5.22) 



Proof: As in the proof of Proposition I5.H we start by locally carrying out the normal 
form transformation in the integral defining the Dirichlet form. Next we define a slightly 
different neighbourhood of the saddle, 



j=q+l 



(5.23) 



where we now choose 



, j = q + l,...,d , 



(5.24) 



with 6 = v^Ke|loge|. The reason for this choice is that we want the potential to be smaller 
than —6"^ on the "sides" {±^1} x of the box. Indeed, we have for all y± G C^, 



2e{d - q) 



€ -K\\oRe\ 



(5.25) 



for sufficiently small e. As a consequence, if h* = hA,B denotes the equilibrium potential, 
Lemma 15.31 and the a priori bound (j5.17p yield 



h%5i,y^) = 0[e 



Clel-«/2-(^+l/2)d[l + i?l(e)]+i?2(e) 

o((e|loge|)(''-^)/('^+^)) i_ 



(5.26) 



e 



_ft:/2+l/4-d/8 



q/2-{c+l/2)d 



[l + |loge| 



while 



h''{-h,y±) = i-0{e^l^) . 

We can now proceed to deriving the lower bound. Observe that 

cap^(5) = ^(^AyjByih*) ^ ^cS^'') ■ 
Now we can write, for any h G TCab, 



(5.27) 
(5.28) 

(5.29) 



ct J -Si 



-V(v)/e( dh{yi,y_L] 
V dyi 



dyi dyj 
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and thus 



(h") >e [ inf e-^(^)/^ f'{yif dyi 

Jet L/: /(±5i)=/i*(±5i,yx) ^-<5i 

The Euler-Lagrange equation for the variational problem is 



dyx . (5.30) 



dV 

ef"{yi)-g^^{yi,y±)f'{yi) = (5.31) 
with boundary conditions h*{—5i,y±) in —6i and h*{5i,y±) in 5i, and has the solution 



fivi) = h*iSi,y±) - [h'iSi,y±) - h^i-Si,y±)]^^ . (5.32) 



<5i 



As a consequence, 



P^(3/i.f-L)/e 

/'(yi) = [h^i6i,y^) - h*{-5uy±)]—^ , (5.33) 



so that substitution in (j5.30|) yields 



-Si 



l2 



^c,(" ) ^ ^ 1 ^y-^ ■ (^•^^) 



The bounds (j5.26p and (j5.28p on h* show that the numerator is of the form 1 — 0(e^/^). 
It now suffices to use the normal form ()5.ip of the potential, and to carry out integration 
with respect to yq+i, ■ ■ - yd- □ 

5.3 Non-quadratic saddles 

Proof of Theorem 13. 2[ For the upper bound, it suffices to apply Proposition 15.11 in 
the case r = 4, g = 2, lii(yi) = C^yf and M2(y2) = A2y|/2. The conditions for the upper 
bound are fulfilled for = (de|loge|/C4)^/^ 62 = 2((ie|loge|/A2)^/2 ^nd c = 0. This yields 
error terms Ri{e) = 0{e^^'^\loge\^^'^) and R2{s) = C'(e'^/^+-^/|loge|^/^). The integrals over 
yi and 7/2 can be computed explicitly as extending their bounds to ±00 only produces a 
negligible error, and we see that the contribution of R2{£) is also negligible. A matching 
lower bound is obtained in a completely analogous way, using Proposition 15.41 with (5i, 
82 of the same order as 5i, 82-, respectively, yielding that R-i{e) is of the same order as 
Ri{e). □ 

Proof of Theorem 13.51 The proof is essentially the same as the previous one, only 
with the roles of 81 and 82 interchanged. □ 
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Proof of Theorem 13. 7[ We again apply Propositions [O and 15.41 now with g = 3, 
ui{yi) = |Ai|yf/2 and ^2(2/2,2/3) = V4:{y2,y3)- The conditions for the upper bound are 
fulfilled for 61 = (2(ie|loge|/|Ai|)i/2, S2 = {2ds\loge\/K_y/'^ and c = 0, and similarly 
for the lower bound. This yields error terms Ri{e) = ©(e^/^llog ej^/^) and R2{s) = 
C(e('^+i)/V|loge|). The integral over yi can be computed explicitly. Writing the integral 
over 7/2 and 7/3 in polar coordinates, and performing the integration over r yields the stated 
expression. □ 



5.4 Bifurcations 

We decompose the proof of Theorem 14.21 into several steps, dealing with positive and 
negative A2 separately. 

Proposition 5.5. Under the assumptions of Theorem\4.2\ and for A2 > 0, 



-P.(B,^^^MM.^/-'/.e-'WA[l. . (5.5, 

where i?+(e) is defined in (|4.13|) . Ia,e is the integral 

/•oo 

Ia,e= / e-(^'+"^')/2.^^^ (5^3g^ 



00 



and a = A2/\/2C4. 

Proof: It suffices to apply Propositions 15.1] and [5741 with r = 4, q = 2, ui a quadratic 
function of yi and U2 a polynomial of degree 4 in y2 and c = 0. We only need to take 
some care in the choice of the 6i. The conditions yield 61 = (2(ie|loge|/|Ai|)^/^ and 



2 _ -A2 + VA| + 32dC4e|loge| 



4G 



4 



(5.37) 



For A2 > (e|loge|)^/^, this implies that 62 is of order (e|log e|/A2)^/^, while for < A2 ^ 
(e|log el)"*^/^, it yields 82 of order (e|log e|)^/''. The expressions of 5i and 82 are similar. 
This yields the stated error terms. The integral over yi is carried out explicitly, while the 
integral over y2 equals /a^e/(2C4)^/^, up to a negligible error term. □ 

Proposition 5.6. Under the assumptions of Theorem\4-^ and for A2 < 0, 



where i?„(e) is defined in (|4.13|) . J^^^ is the integral 

J,,= re-(^^-V4)V2^dx, (5.39) 



and b = ^2/\/2C4. 
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Proof: First note that for small negative A2, 

2 

U2{y2) = Cil 



yl 



8C4 



+ 0{\\2\"^yl) 



(5.40) 



where the constant term corresponds to V{z±) — V{z) (recall z = 0). The situation is more 
difficult than before, because U2 is not increasing on R_|_. When applying Proposition 15. H 
we distinguish two regimes. 

• For /i2 < (e|log e|)^/^, it is sufficient to choose 62 of order (e|loge|)-^/^. 

• For /i2 ^ (e|loge|)^/^, we cannot apply Proposition 15.11 as is, but first split the inte- 
gral over 7/2 into the integrals over ]R+ and over M_. Each integral is in fact dom- 
inated by the integral over an interval of order (e|loge|//X2)^^^ around the minimum 
2/2 = i(/^2/8C4)^/^, so that one can choose 82 of that order. 

We make a similar distinction between regimes when choosing 82 in order to apply Propo- 
sition [531 This yields the stated error terms, and the integrals are treated as before. □ 



Proof of Theorem 14. 2[ In order to complete the proof of Theorem 14.21 it remains to 
examine the integrals Ia,e and Jb^e- First note that 



'27reV2 
1 + a 



where a = aj^fe and 



1 + a 
27r 



,-(S/4+QS/2)/2 



dy 



The change of variables y = z/ vT~-Fa yields 



expi -- 



+ 



az 



1 + a 



dz 



(5.41) 



(5.42) 



(5.43) 



which allows to show that ^'-i- is bounded above and below by positive constants, and to 
compute the limits as a ^ and a ^ 00. The expressions in terms of Bessel functions 
are obtained by observing that 



exp 



+ 28y^ + 



dy 




(5.44) 



because it satisfies the equation f"{8) = {8'^/4)f{8). The other integral is treated in a 
similar way. 



□ 



Proof of Theorem 14. 4[ For Ai < 0, the proof is analogous to the proof of Theorem l4.21 
with the roles of yi and 1/2 interchanged. The same applies for positive Ai up to order 
y/e\loge\. 

For larger Ai, Propositions 15.51 and 15.61 have to be slightly adapted: 

• For the upper bound, we define a neighbourhood of z+ and a neighbourhood 
of Z- in the usual way. Instead of two regions D± , we construct three regions D_ , Dq 
and D+, intersecting respectively 0-,0o and 0+ and contained in the corresponding 
basins of attraction (see Figure [T^) . They are separated by layers Sf^ . The function 
/i+ is then defined to be equal to 1 in D^, to 1/2 in Dq and to in D^. Inside the 
boxes C^, /i+ is constructed in a similar way as before, only with different boundary 
conditions. This yields a factor 1/2 in the capacity. 
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• For the lower bound, we first construct boxes around the saddles in the same 
way as before. Then we connect C'^ and C~ by a tube staying inside Oq (whose 
cross-section is of the same size as the sides of the boxes). One can define coordinates 
(yi, . . . , ?/rf), given by the normal-form transformation inside the boxes, and such that 
yi runs along the length of the tube, in such a way that the Jacobian of the transfor- 
mation X I— > y is close to 1. Lemma 15.31 is still applicable, and yields a priori bounds 
on the equilibrium potential on the sides of the boxes not touching the tube (that is, 
contained in 0±). The Dirichlet form is then bounded below by restricting the domain 
of integration to the union of the boxes and the connecting tube. The remainder of 
the proof is similar, except that the range of yi is larger. The value of the integral is 
dominated by the contributions of the two boxes. □ 

Proof of Theorem 14. 7[ We first consider the case A2 ^ 0. It suffices to apply Propo- 
sitions 15.11 and 15.41 taking some care in the choice of the 5i. The conditions yield 

5i = (2(ie|loge|/|Ai|)^/2 

2 -A2 + \/A| + 32dJC„£|log£| 

^2 = • (5.45) 

For A2 > (e|loge|)^/^, this implies that 82 has order (e|log e|/A2)"'^/^, while for ^ A2 ^ 
(e|log el)"*^/^, it yields 82 of order (e|log el)"*^/^. The expressions of 5i and ^2 are similar. 
This yields the stated error terms. The integral over yi is carried out explicitly, while the 
integral over 2/2 and gives, using polar coordinates. 



U2{y2,yz)/e dy2dy3= r'e-(^2r2+2fc(^)r4)/2e^ ^^^^ 



JBsJO) Jo Jo 



= / \ e-(^'+-(^)^'^)/2ydyd^, (5.46) 

Jo V^k{ip) Jo 

where a{(f) = X2/ \/2ek{(p). Using the fact that 

f e-(^^+'^)V2 ydy = l e'^^/^ dz = ^^-d) , (5.47) 



a straightforward computation shows that the integral over y is approximated by 

e+(a(v?)) 

l + a{ip) 

For small negative A2, we can write 

I \ I \ 2 \2 

2 IA2I \ Ar) 



(5.48) 



where the constant term corresponds to the actual minimum of the potential. 

If — (e|loge|)^/2 < A2 < 0, it suffices to apply Propositions 15. II and 15.41 with 82 of order 

{e\loge\)'/\ □ 

Proof of (|4.42p . The main task is to compute the integral related to U2, namely 

J :=l e-["^^'+2Er^2 C^.r^'V^e jT'" e-I)2Mr-^ cos(2M^)/. . (5.50) 
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We first carry out the integral over ip, yielding 



^ ^ ^-D^^^r-M ,os(2M^)/. ^ 27r/o (^^^^) > (5-51) 



where Iq is the modified Bessel function 



h{a) = ^ jJe^'^'^P dip. (5.52) 

The Laplace method shows that for large a, /o(ct) behaves like e" /\/27ra. We thus intro- 
duce the bounded function 



vr 



27r 



Inserting in (j5.50p and performing the change of variable = (e/2C4)^/^n yields 

J = ^ f ^~[u^-2{\\2\/Vm)u+0{u^)]/2e X{D2Me~^{u/{2Cifl^f^) 

V8C-J Vl + ^2A/e-i(V(2C4)i/2)M 

Applying the Laplace method shows that the integral is dominated by u close to u* = 
|A2|/(8C4)i/2. Relating the obtained expression to the eigenvalues at the new saddles via 
the relation 

D2M ( U* \M _ fi2f^3 1 + 0{X2) , . 

e \{2CiY'^) (2M)2 8eC74 ^^'^^^ 
yields the necessary control on J . The stated formula for the capacity follows. □ 



A Normal forms 

Proof of Proposition 12.81 Let us denote by Qk{n,m) the vector space of functions 
g: M" — > M"* which are homogeneous of degree k (i.e., g{tx) = t^g{x) Vx). We write the 
Taylor series of V in the form 

V{x) = V2{x) + y3(x) + Vi{x) + o{\\x\\t) , (A.l) 

where Vk G Qk{d,l) for k = 2,3,4. We first look for a function (72 G Q2{d,d) such that 
V o [id +(72] contains as few terms of order 3 as possible. The Taylor series of y o [id +52] 
can be written 



V{x + g2{x)) 

= V2{x) + VV2{x) ■ g2{x) + V^{x) + V2{g2{x)) + VV^{x) ■ g2{x) + Vi{x) +o{\\x\\t) . 

V ' ^ V ' 

order 3 order 4 

(A.2) 

Now consider the so-called adjoint map T: Q2id, d) — > ^3(d, 1), (72 ^ VV2(-) • 52, seen as a 
linear map between vector spaces. All terms of V^{x) in the image of T can be eliminated 
by a suitable choice of 52- Let ei denote the /th vector in the canonical basis of IR*^. We 
see that 

T{xjXkei) = XiXjXkXi / for / = 2, . . . , d . (A. 3) 
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Thus all monomials except x'^ are in the image of T. Since T involves multiplication by 
X2 or X3 or . . . or Xd, however, xf is not in the image of T. Hence this term is resonant. 
We can thus choose g2 in such a way that 

V{x + g2{x)) = V2{x) +ynixl + y2{g2{x)) + VV^jx) • g2{x) + V4{x) +oi\\x\\l) . (A.4) 

order 3 order 4 

Now a completely analogous argument shows that we can construct a function g^ S (d, d) 
such that V o [id +52] ° [id +53] has some constant times xf as the only term of order 4. It 
remains to determine this constant. From (jA.4p we deduce that it has the expression 



-, d d 

i=2 j=i 



(A.5) 



where g^^ denotes the coefficient of xfe^ in g2- The expression of T shows that necessarily 
-Viij/Xj for j = 2, . . . , d, while we may choose gl^ = 0. This yields ()2.14p . □ 
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